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ABSTRACT 


An  analysis  of  fin-line  characteristics  is  presented. 

The  dispersion  characteristic  and  characteristic  impedance 
of  the  structure  are  determined  by  the  application  of  the 
spectral  domain  technique.  Numerical  results  are  presented 
and  compared  with  known  data  for  ridged  waveguide,  dielectric 
slab  loaded  waveguide  and  slot line.  Excellent  agreement  has 
been  obtained.  Design  curves  are  also  presented  for  practi¬ 
cal  millimeter  wave  fin-line  structures  covering  the  26.5  - 
140  GHz  frequency  range. 
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INTRODUCTION 


I . 


A.  BACKGROUND 

A  variety  of  waveguiding  structures  have  been  proposed 
for  use  at  microwave  and  millimeter  wave  frequencies  [Ref.  1] . 
In  recent  years,  one  of  these  structures,  fin-line,  has  gained 
in  importance  as  a  transmission  medium  in  millimeter  wave 
circuit  construction  [Refs.  2-7].  Fin-line  has  been  found 
superior  to  microstrip  at  millimeter  wavelengths  as  the 
former • provides  eased  production  tolerances,  better  compati¬ 
bility  with  hybrid  devices,  greater  freedom  from  radiation 
and  higher  mode  propagation,  combined  with  the  ability  to 
construct  simple  transitions  to  conventional  rectangular 
waveguide. 

Figure  1  shows  a  3-dimensional  view  of  fin-line.  This 
structure  is  interesting  in  that  it  may  be  viewed  in  various 
ways  depending  upon  the  value  of  W/b,  the  ratio  of  slot  width 
to  waveguide  height.  For  small  values  of  W/b,  the  structure 
may  be  appropriately  viewed  as  a  slotline  with  a  rectangular 
shield.  For  large  values  of  W/b,  the  structure  is  more  easily 
seen  as  a  ridged  waveguide  loaded  with  dielectric.  When 
W/b  =  1,  a  dielectric  slab  loaded  waveguide  exists  and  finally 
for  W/b  =  1  and  er  =  1,  the  structure  takes  on  the  configura¬ 
tion  of  empty  rectangular  waveguide. 

The  most  commonly  used  transmission  line  to  date  in 
microwave  integrated  circuits  is  the  microstrip  line.  This 
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structure  when  used  at  high  frequencies  (millimeter  wave 
range)  results  in  problems  which  include  radiation  loss, 
spurious  coupling,  dispersion  and  higher  mode  propagation. 

In  fin-line  structures,  metal  fins  are  printed  on  a 
dielectric  substrate  which  bridges  the  broad  walls  of  a 
rectangular  waveguide.  In  effect,  the  line  is  a  printed 
ridged  waveguide  and  can  be  designed  to  have  a  wider  useful 
bandwidth  than  conventional  line  and  can  provide  bandwidths 
in  excess  of  an  octave  with  less  attenuation  than  microstrip. 
This  adaption  of  ridged  waveguide  permits  circuit  elements 
to  be  fabricated  by  photo-etching  at  low  cost  and  is  compati¬ 
ble  with  thin-film  hybrid  techniques. 

In  passive  circuits,  such  as  filters,  the  fins  as  illus¬ 
trated  in  Figure  1  may  be  directly  grounded  to  the  metal  walls 
of  the  waveguide,  and  lumped  elements,  such  as  beam- lead 
capacitors,  may  be  added.  The  gap  between  the  fins  can  be 
varied  along  the  longitudinal  axis  to  provide  low-cost  circuit 
elements.  When  fin-lines  are  constructed  as  illustrated  by 
Figure  2,  the  waveguide  is  parted  along  a  plane  where  the 
current  flow  is  parallel  to  the  break,  as  in  a  common  slotted 
line  and  the  dielectric  material  extends  through  the  separation 
in  the  broad  wall  of  the  shield.  The  upper  fin  is  insulated 
from  the  housing  at  dc  by  a  dielectric  gasket.  By  making  the 
thickness  of  the  broad  walls  a  quarter  wavelength,  a  short 
circuit  will  appear  between  the  fin  and  the  inner  wall  of  the 
shield  at  RF.  Thus  a  practical  fin-line  will  have  the  same 
electrical  characteristic  as  the  idealized  structure  of 
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Figure  1.  The  lower  fin  is  grounded  directly  by  a  metal 
gasket  to  provide  a  dc  return  when  solid  state  devices  are 
mounted  on  the  fins. 

B .  RELATED  WORK 

Fin-line  was  first  proposed  in  1972  by  Meier  [Refs.  2 
and  3]  and  has  since  received  considerable  attention  by  a 
number  of  other  investigators.  Meier's  work  was  primarily 
experimental  in  nature  and  the  wavelength  and  impedance 
calculations  were  approximated.  Previous  analytical  work  on 
fin-line  has  been  based  on  various  methods  [Refs.  8-11]. 

In  this  study,  dispersion  of  different  fin-line  geometries 
was  obtained  by  applying  Galerkin's  method  to  the  electric 
fields  in  the  slot  interface  [Ref.  8] .  Another  approach  used 
the  Transmission  Line  Matrix  (TLM)  technique  yielding  the 
dispersion  characteristic  of  the  fundamental  and  higher  order 
modes  of  propagation  [Ref.  9]. 

C.  STATEMENT  OF  OBJECTIVES 

The  Spectral  Domain  technique  is  used  here  in  the  analy¬ 
sis  of  the  fin-line  structure  shown  in  Figure  1.  This  tech- 
nicue  was  first  suggested  by  Itoh  and  Mittra  [Ref.  12]  and 
has  since  been  applied  to  the  analysis  of  other  structures 
[Refs.  13-15]. 

One  of  the  advantages  of  this  approach  is  that  it  is 
numerically  more  efficient  than  the  conventional  methods  that 
work  directly  in  the  space  domain.  This  is  due  primarily  to 
the  fact  that  the  process  of  Fourier  transformation  of  the 
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coupled  integral  equations  in  the  space  domain  yields  a  pair 
of  algebraic  equations  in  the  transform  domain  that  are 
relatively  easier  to  handle.  Another  important  advantage  is 
that  the  Green's  function  takes  a  much  simpler  form  in  the 
transform  domain,  as  compared  to  the  space  domain  where  no 
convenient  form  of  the  Green's  function  is  known  to  exist. 

The  Spectral  Domain  Transform  method  is  first  applied 
here  in  the  formulation  of  the  dispersion  characteristic  of 
the  fin-line  structure  of  Figure  1.  During  the  application, 
a  new  matrix  approach  to  the  implementation  of  this  method 
is  presented.  The  characteristic  impedance  is  next  formu¬ 
lated  using  this  technique  in  terms  of  the  dispersion  char¬ 
acteristic.  Following  the  theroetical  analysis,  an  explana¬ 
tion  of  the  computer  program  used  in  determination  of  the 
wavelength  and  characteristic  impedance  is  presented.  Numeri¬ 
cal  results  are  then  compared  with  known  data  for  ridged 
waveguide,  slab  loaded  waveguide  and  slot line.  These  com¬ 
parisons  establish  the  accuracy  of  the  numerical  results  and 
illustrate  the  applicability  of  this  method  for  the  full 
range  of  structure  parameters.  Finally,  several  families 
of  design  curves  are  presented  for  a  practical  choice  of 
fin- line  parameters. 
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II.  THEORETICAL  ANALYSIS  OF  FIN-LINE 


A.  FIELD  AND  BOUNDARY  CONDITIONS 

Fin-line  supports  a  hybrid  field  and  it  is  known  that  all 

hybrid  field  components  can  be  obtained  from  the  superposition 

of  TE  and  TM  modes  which  are  related  to  the  two  scalar  poten- 
0  h 

tial  functions  4>  (x,y)  and  $  (x,y)  ,  where  the  superscripts  e 
and  h  denote  electric  and  magnetic,  respectively  [Ref.  16]. 

The  axial  components  of  TM  and  TE  are  then 

E  =  k l  <(>e (x,y)  erz  (1) 

C 

h2  =  kc  *h(x'y}  eFZ  (2) 

2  2  2  2 

where  k  .  =  k.  +  r  and  k.  =  w  u.e.,  i  =  1,  2,  3  for  each 

CX  X  X  XX 

of  the  three  regions  defined  in  Figure  (1). 

Through  Maxwell's  curl  equations  the  transverse  field 
components  are  then  determined  by  these  axial  components  and 
can  be  given  as 
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(6) 


Since  the  lossless  case  is  assumed,  the  propagation  constant 

2  2  2 

is  r  =  ±je  and  k  .  =  k.  -  8  .  We  will  also  assume  that 

Cl  X 

e,  =£..  =  £  and  £_=££. 

13o  2  o  r 

Applying  boundary  conditions  at  the  walls  in  region  1, 
tangential  field  components  must  be  zero,  so  it  follows  that 
at  y  =  h^  +  D 

3H  . 

Ezl(x,hi+D,z>  =  0  ^i(x,h1+D,z)  -  0  (7) 


Exl(x,hl+D,2)  =  °  -3^(x,hi+D,z)  =  0.  (8) 

At  the  interface  between  region  1  and  2,  tangential  field 
components  must  be  continuous. 

At  y  =  D 


E  ,  (x, D,z)  =  E  - (x,D, z)  (9) 

zl  z2 

Exl(c,E),z)  =  Ex2<x,D,z)  .  (10) 

Also  the  electric  fields  at  y  *  D  will  exist  only  in  the 
slot  and  can  be  expressed  as, 
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Hz2 (x,0,z) 


Hz3(x,0,z) 


(17) 


Hx2(X,0'Z)  =  hx3(x»°'z)  • 


(18) 


Once  again  at  the  shield  wall  in  region  3,  the  tangential 
field  components  must  be  zero,  therefore  at  y  =  -h_ 


3H 


Ez3(x,-h2,z)  =  0 


Ex3(x,-h2,z)  =  0 


z3 


3y  (x, ~h2 / z)  =  0 


3H 


x3 


=  0 


(19) 


(20) 


The  final  boundary  conditions  occur  at  x  =  ±b/2  where  the 
tangential  components  must  be  zero  in  all  regions.  At 
x  =  ±b/2 


Ezi(±b/2,y,z)  =  0 


Exi(lb/2,y,2)  *  0 


(21) 


(22) 


B.  SPECTRAL  DOMAIN  APPROACH  TO  DISPERSION  CHARACTERISTIC 
The  approach  used  in  the  analysis  of  the  dispersion 
characteristic  is  basically  a  modification  of  Galerkin's 
approach  adapted  for  application  in  the  Fourier  transform 
or  spectral  domain. 

The  scalar  potential  functions  $e  and  satisfy  the 
Helmholtz  equations  in  the  three  spatial  regions,  thus 
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0 


(23) 


2 

where  v  denotes  the  two-dimensional  Laplacian  operator  in  the 
transverse  direction.  These  equations  are  second-order  differ¬ 
ential  equations  for  the  unknown  potential  functions  which 
may  be  Fourier  transformed  with  respect  to  x  to  obtain  ordinary 
differential  equations  via  the  transformation  defined  as 

“  j  a  X 

F{<j>(x,y)}  =  $(an/y)  =  /  <J>(x,y)  e  n  dx  .  (24) 


The  Helmholtz  equations  of  (23)  are  transformed  into 


3  ^“n'7* 


kci> 


(25) 


0  Vt 

for  the  electric  field,  ,  and  the  magnetic  field,  $  . 

The  solutions  to  these  ordinary  differential  equations 
for  the  three  regions  after  applying  the  boundary  conditions 
at  y  =  h^  +  D  and  y  =  -h2  (equations  (7),  (8),  (19)  and  (20)) 
can  be  written  as 


(an»y)  *  Ae(an)  sinh  Y1(D+h1~y)  (26a) 

«|(an»y)  =  Be(an)  sinh  Y2y+Ce(an)  cosh  y  2y  (26b) 
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(26c) 


*3(Vy) 
^l *an,y^ 
*2(an'y) 
*3(an'y) 


De(a  )  sinh  Y,(h*+y) 
n  3  2 

Ah(an)  cosh  y^dHh^-y)  (26d) 

Bh(a  )  sinh  Y-,y  +  C^(a  )  cosh  Y-,y  (26e) 
n  2  n  2 

Dh(an)  cosh  Y3(h2+y)  (2 6f) 


where 


(27) 


and 


n2?r/b  even 

(2n-l)  7r/b  $h  odd 


(28a) 


(28b) 


.  .  2 
It  is  important  to  observe  at  this  point  that  y^  may  be 

less  than  zero  in  any  of  the  three  regions  of  the  structure 

under  certain  conditions.  Therefore  three  sets  of  solutions 

can  exist  for  equations  (26a)  to  (26f) .  When  afl  =  0  and 

approaches  kQ  (where  kQ  is  the  wave  number  for  free  space) 

2 

8  is  less  than  k^  and  so  y.  <  0.  Under  this  condition  the 

hyperbolic  functions  in  all  three  regions  are  replaced  by 

.  .  2  2 
trigonometric  functions.  If  kQ  <  s  <  k2,  then  and  y3 
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2 

are  greater  than  zero  and  <  0  for  some  values  of  ^  and 

the  trigonometric  functions  replace  the  hyperbolic  functions 

in  the  spatial  region  2  only.  This  suggests  that  the  nature 

of  the  field  is  dependent  upon  the  values  of  the  transform 

2 

variable,  a  .  The  last  solution  occurs  when  y.  >  0  for  which 
n  i 

2 

the  solutions  are  given  above.  For  conditions  when  y^  <  0, 
y ^  replaces  y^  such  that  (y^j  =  -yf  and  cosh  y^y  is  replaced 

H  It 

by  cos  Yiy  and  sinh  y^y  is  replaced  by  j  sin  y^y. 

The  next  step  in  obtaining  a  solution  to  the  differential 

equations  is  the  application  of  the  continuity  or  boundary 

conditions  at  the  interfaces  between  the  three  spatial  regions, 

namely  at  y  =  0  and  y  =  D.  These  conditions  are  given  in 

equations  (9)  through  (18)  .  Through  their  application,  the 

0  h 

eight  coefficients  A  through  D  may  be  related  to  each 
other,  to  the  fin  surface  current  J(x)  and  the  slot  field 
e(x).  The  resulting  set  of  linear  equations  (really  a  set  of 
transformed  boundary  equations)  may  be  written  in  a  matrix  form 
as  follows: 

0 
0 

• 

Ex 
Ez 
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where  E. ( a  )  is  the  transform  of  the  slot  field  and  J. (a  ) 
in  i  n 

is  the  transform  of  the  surface  current.  The  matrices  [Mg] 

and  [Mj]  are  square  8x8  matrices  and  differ  only  in  the 

last  two  rows.  Details  of  the  application  of  the  continuity 

conditions  and  the  formation  of  the  [Mg]  and  [Mj]  matrices 

are  shown  in  Appendix  A,  with  a  listing  of  the  matrix  elements 

2  2 

under  the  conditions  y^  >  0  and  y^  <  0  and  normalized  with 
respect  to  D. 

To  solve  the  matrix  equations  (29) ,  the  unknown  coeffi- 
•  g  h 

cients  A  through  D  are  determined  by 


where  [Mg]  ^  is  the  inverse  matrix  of  [Mg] .  Now  substituting 


eauation  (30)  into  (29b)  results  in  the  single  matrix 
equation 


0 

0 

0 

0 

IMj] [Me]_1 

• 

= 

• 

Ex 

Jx 

.  Jz 

Using  the  four  lower  right  hand  corner  elements  of  the 


Er\  HD 

is  converted 

to 

G.  ( a  ,  8 ) 

I  n 

G2(an'3) 

E  (a  )  1 
x  n 

J  (a  ) 
x  n 

.G, (a  , 8) 
j  n 

G4(an,3) 

E  (a  ) 
z  n 

J  (a  ) 

z  n 

(32) 


where  the  elements  of  the  matrix  [G]  are  the  Fourier  transforms 
of  the  components  of  the  dyadic  Green's  function  for  this 
structure . 

A  solution  to  equation  (32)  is  obtained  by  using  the 
Method  of  Moments  [Ref.  17] .  The  inner  product  is  defined  as 


00 

<f(ctn),  g(an)>  =  l  f(o  )  g*Un)  .  (33) 

n=-«  11 
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Galerkin’s  method  is  used  here  and  thus  the  Fourier  trans¬ 
forms  for  the  electric  field  components,  E^  and  E will  be 
used  for  the  weighting  functions,  (1/(1^).  Taking  the  inner 
product  of  equations  (32)  with  ®  an<* 

U)  (a  )  =  E  (a  )  yields 
2.  n  z  n 

<GlEx'Ex>+<G2Ez'Ex>  =  ° 

(34) 

<G3Ex'E2>+<G4Ex'Ez>  "  0  * 

Following  from  Parseval's  theorem,  the  right-hand  side  of 
these  equations  are  zero  because  of  the  orthogonality  of 
the  electric  field  and  surface  current  at  y  =  D. 

Up  to  this  point,  the  formulation  of  the  problem  is  exact 
since  no  approximations  have  been  made.  The  electric  field 
components  and  their  transforms  although  unknown  can  be 
expanded  in  a  set  of  basis  functions.  An  investigation  of 
various  one-term  approximations  made  by  hnorr  and  Kuchler 
[Ref.  18]  suggests  good  accuracy  can  be  achieved  if  the 
electric  field  distribution  between  the  fins  is  approximated 
by 


ex(x) 


1  I  x  I  <_  W/2 

0  elsewhere 


ez(x)  =  0 


(35) 
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for  the  dominant  mode.  This  choice  (shown  in  Figure  3  for 
A  =  1)  has  been  shown  to  give  accurate  results  for  slotlines 
with  W/D  <  2  (Ref.  18]  and  it  is  exact  for  dielectric  slab- 


loaded  waveguide  (W/b  =  1) .  Assuming  that  e  (x)  =0  reduces 

z 

equation  (34)  to 


<GiEx'Ex>  -  l  ‘Vv^iW'2  *  0  1361 

n=-°° 

where  E  (a  )  is  the  transform  of  the  electric  field  distri- 
x  n 

bution 

W/2  ja  x  sin (a  W/2) 

Ex(an>  *  _M//  ex(x>  e  "  dx  =  AW  an  W/2 -  . 

(37) 

So,  one  can  see  that  (36)  is  dependent  on  the  ratio  of  X/X' 

or  effective  dielectric  constant,  e  ,  defined  by 

ref  f 

e  =  (X/X')2  =  62/k2  . 

ref  f 

The  dispersion  characteristic  of  the  dominant  fin-line 
mode  can  now  be  found  by  varying  6/  the  propagation  constant, 
such  that  equation  (36)  is  satisfied  for  a  given  set  of 
physical  parameters  at  a  desired  frequency  of  operation. 


4 


-w/2 


Figure  3 


♦  w/2 


+  X 


Assumed  Electric  Field  Component  in  Slot 
in  x-direction  Versus  x  for  Fin-line 
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C.  CHARACTERISTIC  IMPEDANCE  IN  TERMS  OF  DISPERSION 

CHARACTERISTICS 

The  definition  of  the  characteristic  impedance  for  an 
ideal  TEM  transmission  line  is  uniquely  given  by  static 
quantities.  Since  fin-line  is  not  a  TEM  transmission  line 
nor  does  it  support  a  pure  TE  or  TM  waveguide  mode  but  rather 
supports  hybrid  modes,  no  unique  definition  of  the  character¬ 
istic  impedance  can  be  found. 

One  possible  choice  is  to  define  it  as 


z 


o 


o 

2  P 

avg 


(33) 


where  Vq  is  the  slot  voltage  defined  as 


V 


/E  •  dz  . 

slot 


Assuming  the  electric  field  shown  in  Figure  3  and  (35) , 


W/2 

V  =  f  A  dx  =  1 

°  -W/2 


where  A  was  arbitrarily  selected  as  1/W  so  that  AW  =  1. 

Pavg  *s  t*le  time-averaged  power  flow  in  the  fin-line  struc¬ 
ture  which  is  given  by 
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(39) 


P  =  4  Re  /  /  E  x  h*-  a  da 

avg  2  1  J  2 

s 

=  \  Re  /  /  (E  H  *  -  E  H  *)  dx  dy  . 

2  J  }  x  y  y  x 
s 

Substituting  equations  (3) ,  (4) ,  (5)  and  (6)  into  (39) 
results  in 


•  k  Re  / 


-b/2 


-  (jB|^-+jUu|^-)(-jB||--ja.e|^-)]dxdy.  (40) 


This  formulation  applies  to  the  time-averaged  power  flow 
in  the  spatial  domain.  Through  the  use  of  Parseval's  theorem, 
it  is  possible  to  transform  this  expression  into  the  spectral 
domain  where  the  scalar  potential  functions  are  known.  The 
integration  with  respect  to  x  is  transformed  via  Parseval's 
relation  into  an  infinite  summation.  Applying  this  theorem 

0  H 

where  *  (x,y)  and  $  (x,y)  are  real  yields 
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Therefore,  equation  (40)  is  transformed  into 


^  I  a 

1  D  r°°  (  0  2 1  e  1 2  2 1  h ,  2  ,  3$e 

2b  Re  l  I  {-Buea  |4  |  -8«ua  |  4  |  -Sa)e|  — 

n=-oo  -h^  “ 1 


,  3  <th  i  2  .  „2  fie  3(4h)*  ,  3  4e  (1h, 

-  Bow  I  ay-1  -  36  «„(♦  gy  +  3y"  ( *  >*] 


g  h 

,  ..2  fAh  3(4  )*  ,  3  4  ,  .  e.  *  •,  i  , 

+  D|c  an[$  S—L.  +  —(4  )*]}  dy. 


This  expression  must  be  evaluated  in  each  of  the  three 
regions  of  the  fin-line  shown  in  Figure  1.  Therefore  the 
power  flow  may  be  expressed  by 


P  *  P.  +  P-  +  P, 

avg  12  3 


where  P^  P2  and  P3  are  defined  as 


55  Re  l  /  {f(4®,4j)>  dy 

n=-°=  d 


(42a) 


00  Q 

-  -k  Re  i  /  }dy 


n=-»  o 


(42b) 


^  Re  l  /  {f(4®,45)>dy 

"  n=-«  -hj 


(42c) 
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The  simple  functional  dependence  on  the  variable  y  allows 
one  to  perform  the  integration  in  equations  (42a) -(42c) 
analytically  (see  Appendix  B) .  This  leaves  Pavg  as  the 
summation  of  the  power  flow  in  each  region. 

Pavg  *  1  <P1  +  p2  +  P3'  '  <431 

n=-°° 

Now  that  ?aVg  can  be  determined  by  the  application  of  the 
dispersion  characteristic  or  propagation  constant  determined 
earlier  and  that  the  slot  voltage  Vq  is  known,  the  character¬ 
istic  impedance  can  be  computed.  The  details  of  these  lengthy 
but  straightforward  algebraic  manipulations  are  shown  in 
Appendix  B. 


III.  COMPUTER  PROGRAMMING 

In  Section  II  the  theoretical  background  of  the  analysis 
was  provided.  This  section  concerns  itself  with  the  computer 
programming  technique  used  for  fin-line  analysis.  The  program 
was  designed  for  the  determination  of  the  wavelength  and  the 
characteristic  impedance  of  a  particular  fin-line  structure 
over  a  range  of  frequencies  but  due  to  the  nature  of  the 
structure  various  other  configurations  may  be  considered  by 
the  proper  selection  of  input  parameters.  Dielectric  slab 
loaded  and  ridged  waveguide  along  with  rectangular  waveguide 
may  also  be  analyzed.  Before  che  program  is  described  some 
general  considerations  which  lead  to  its  development  are 
necessary  for  understanding  of  the  numerical  methods  used. 

A.  REMARKS  ON  NUMERICAL  ANALYSIS 

The  computation  of  the  characteristic  impedance  is  based 
upon  the  solution  to  the  dispersion  characteristic  problem 
for  the  transmission  line  under  consideration.  In  other 
words,  the  wavelength  ratio,  \'/\,  or  the  wave  propagation 
constant,  8,  must  be  known  before  any  other  investigations  can 
be  started  since  only  in  this  case  are  the  scalar  potential 
functions  in  the  transform  domain  known.  Hence,  in  the  analy¬ 
sis  of  the  various  structures,  the  starting  point  is  a  solu¬ 
tion  for  8,  the  propagation  constant  of  equation  (36) . 
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In  general  a  8  is  sought  such  that 


l  G,  (a  ,B)  |E(a  )  |2  =  0  (44) 

L  in  n 

n=-» 

where  G^(an,6)  is  a  component  of  the  transformed  dyadic 
Green's  function  and  E(an)  is  the  approximation  to  the 
Fourier  transform  of  the  electric  field  in  the  slot  in  the 
plane  y  =  D.  A  first  order  approximation  is  assumed  here 
which  eliminated  the  other  terms  of  equation  (34) . 

It  is  obvious  that  an  analytic  solution  for  s  cannot  be 
obtained  due  to  the  algebraic  complexity  leading  to  the  quan¬ 
tity  G^ (°n»  6)  as  shown  in  Appendix  A.  Thus  a  numerical 
solution  is  mandatory  and  must  be  obtained  by  computer 
programming . 

It  is  shown  in  Section  II  that  the  transforms  of  the  dyadic 
Green's  function  components  were  obtained  by  mathematical 
manipulation  of  matrices.  The  matrix  elements  are  either 
real  or  imaginary  so  that  in  general  a  complex  computation  on 
the  computer  is  necessary.  Furthermore,  it  is  shown  that  the 
elements  change  from  hyperbolic  to  trigonometric  functions. 

H 

For  examaple  sinh  y2d  becomes  jsin  y2d-  was  f°und  that 

these  quantities  are  not  only  pure  real  or  pure  imaginary 

but  change  from  real  to  imaginary  or  vice  versa  for  the 
2  2 

cases  y^  >  0  and  y^  <  0. 

Previously,  equation  (32)  was  obtained  by  lengthy  alge¬ 
braic  manipulation  of  the  boundary  equations .  In  the  approach 
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described  here  the  matrices  (29a)  and  (29b)  are  defined 
directly  from  the  boundary  equations  and  equation  (32)  is 
arrived  at  by  numerical  computation.  There  is  some  sacri¬ 
fice  in  numerical  efficiency  in  this  approach  but  the  formu¬ 
lation  of  the  problem  is  straightforward. 

Since  in  equation  (44)  the  absolute  value  squared  of  the 
assumed  electric  field  is  taken,  this  part  of  the  series  will 
be  in  general  a  real  and  even  function  of  aR.  Preliminary 
numerical  investigations  by  Kuchler  [Ref.  19]  of  the  coeffi¬ 
cients  of  these  series  indicated  an  even  distribution  with 

respect  to  the  variable  a  ,  also  exists.  Therefore,  all  terms 

n 

of  the  series  (44)  are  even  functions  which  allows  the  saving 
of  computation  time  by  summing  over  only  half  the  interval. 

A  desired  solution  for  equation  (44)  can  now  be  found  by 
assigning  an  arbitrary  numerical  value  for  the  wavelength 
ratio  and  then  iteratively  changing  this  until  the  series  con¬ 
verges  to  zero  or  more  precisely  a  small  enough  value  so 
that  a  prescribed  accuracy  for  the  wavelength  ratio  is  realized. 

The  next  task  is  the  determination  of  the  time  average 
power  flow.  Here  again  complex  arithmetic  expressions  are 
involved  as  seen  in  Appendix  B. 

Dielectric  substrates  are  commercially  available  for 
various  thicknesses  D  so  that  in  general  the  line  parameters 
should  be  computed  for  each  thickness  and  frequency.  A  more 
generalized  method  of  analysis  is  simply  to  scale  the  various 
geometric  quantities  involved  with  respect  to  this  thickness 
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to  obtain  more  versatile  parameters  for  a  varying  normalized 
frequency  D/A .  The  resulting  design  curves  may  then  be 
applied  to  the  structures  scaled  to  different  operating 
frequencies . 

B.  COMPUTER  PROGRAM  DESCRIPTION  FOR  FIN-LINE  CONFIGURED 

TRANSMISSION  LINE 

A  computer  program  which  provides  a  numerical  analysis  of 
fin-line  is  shown  in  Appendix  C  with  the  main  program  called 
'FIN-LINE'  and  its  accompanying  subroutines  'CALC'  and  "CMTRIN, ' 
all  written  in  the  FORTRAN  IV  programming  language. 

The  program  'FIN-LINE'  starts  by  reading  the  input  data 
which  includes  the  structure  dimensions,  the  dielectric  con¬ 
stant,  the  initial  values  of  the  slot  width  and  frequency 
and  the  iteration  steps.  After  the  input  data  is  read, 
initialization  of  computation  loops  occur  for  varying  slot 
widths  and  normalized  frequencies.  Next,  the  program  determines 
the  type  of  structure  that  exists  and  whether  or  not  the  con¬ 
figuration  is  operating  in  the  dominant  mode.  Operation  at 
higher  order  empty  waveguide  modes  are  allowed  but  program 
termination  occurs  at  the  onset  of  surface  waves. 

The  program  proceeds  with  preliminary  computations  used 
in  the  determination  of  the  dispersion  characteristic.  The 
first  major  computation  of  the  program  that  takes  place  is 
the  finding  of  the  wavelength  ratio  A '/A  which  satisfies 
equation  (36)  .  The  necessary  coefficients  for  determining 
the  series  are  provided  by  the  subroutine  'CALC'.  An  initial 
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is  the  lower 


value  of  X'/X  =  1/ve  is  selected  where  1/ve 

r  r 

limit  of  x'/x,  that  is  when  x'  =  X  ,  this  wavelength  of  the 

a 

dielectric  material.  The  results  from  the  initial  computation 

are  used  to  select  a  new  value  for  x'/x  through  a  bisectional 

zero  finding  method  to  improve  the  results  and  the  procedure 

starts  over  again.  The  computation  is  terminated  when  the 

correct  value  for  X'/X  is  determined  which  satisfies  equation 

(36)  or  when  the  absolute  value  of  the  difference  between  the 

-4 

old  and  the  new  value  for  X'/X  is  less  than  5*10  and  thus 
produces  a  three  digit  accuracy  after  the  decimal  point. 

Using  the  determined  value  of  the  wavelength  ratio,  the 
program  continues  with  the  computation  of  the  characteristic 
impedance.  The  coafficients  for  this  computation  (equation  40) 
are  provided  again  by  the  subroutine  'CALC'. 

After  completion  of  the  dispersion  and  impedance  compu¬ 
tations,  the  program  prints  these  values  and  restarts  with 
new  values  for  the  normalized  frequency  and/or  slot  width. 

Figure  4  presents  a  flow  diagram  for  the  main  program 
'FIN-LINE'.  All  the  variable  names  have  been  chosen  such 
that  the  statements  can  be  readily  compared  with  the  respec¬ 
tive  equations  of  the  appropriate  appendices. 

The  subroutine  'CALC'  provides  the  coefficients  and  matrix 
elements  as  well  as  an  extensive  amount  of  the  computation 
used  in  the  determination  of  both  the  dispersion  and  impedance. 
The  entry  variable  is  X/X '  for  both  characteristics.  A  flow 
diagram  for  the  function  'CALC'  is  shown  in  Figure  5. 


Figure  4d.  Flow  Diagram  for  Main  Program  'FIN-LINE 
(Continued) 


Figure  4e.  Flow  Diagram  for  Main  Program  ’FIN-LINE' 
(Continued 
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Figure  5a.  Flow  Diagram  of  Subroutine  'CALC' 
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Figure  5c.  Flow  Diagram  of  Subroutine  'CALC' 
(Continued) 


Figure  5f.  Flow  Diagram  of  Subroutine  'CALC' 
(Continued) 
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Figure  5g.  Flow  Diagram  of  Subroutine  'CALC' 
(Continued) 


Figure  5h.  Flow  Diagram  of  Subroutine  'CALC' 
(Continued) 


Figure  5i.  Flow  Diagram  of  Subroutine  'CALC' 
(Continued) 


Figure  5 j .  Flow  Diagram  of  Subroutine  'CALC1 
(Continued) 
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Once  the  subroutine  is  entered  from  the  main  program, 
variables  used  throughout  are  computed.  Following  the  initial 
computations,  a  decision  is  made  as  to  whether  hyperbolic  or 
trigonometric  functions  are  to  be  used  in  the  formation  of 

the  matrix  elements  shown  in  Appendix  A.  This  decision,  of 

.  2 
course,  is  made  by  determining  whether  is  less  than  or 

2 

greater  than  zero.  If  y^  >  0,  hyperbolic  functions  apply 
2 

and  if  y ^  <  0,  trigonometric  replace  the  hyperbolic  functions. 
After  the  decision  making  and  the  formation  of  the  [ME]  and 
[MJ]  matrices,  the  [ME]  matrix  is  inverted  and  the  matrix 
elements  of  concern  in  the  [MJ]  and  [ME]  ^  matrices  are  multi¬ 
plied  together  for  the  formation  of  element  Gil  of  the  dyadic 
Green's  functions. 

The  transform,  of  the  assumed  slot  electric  field  distribu¬ 
tion  is  determined.  At  this  point  the  decision,  as  to  whether 
the  dispersion  or  the  impedance  is  being  determined,  is  made. 
If  the  dispersion  computation  is  in  process,  the  subroutine 
computes  the  terms  of  equation  (36)  and  returns  to  the  main 
program  for  a  new  selection  of  A' /A. 

If  the  characteristic  impedance  computation  is  being  per¬ 
formed,  the  program  continues  with  the  determination  of  the 

.  .  6  h 

coefficient  matrix  elements  A  through  D  .  Once  again  the 

decision  for  the  use  of  trigonometric  or  hyperbolic  function 

is  made  followed  by  additional  variable  computations  needed 

in  determination  of  power. 
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2 

The  power  flow  in  each  region  is  determined  for  y^ 
greater  or  less  than  zero,  depending  on  the  conditions,  and 
totaled.  This  procedure  is  followed  until  the  difference 
between  two  consecutive  values  of  the  impedance  is  less  than 
0.1%  at  which  time  return  to  the  main  program  occurs. 

The  subroutine  'CMTRIN'  takes  the  complex  [ME]  matrix 
and  inverts  it  into  the  [ME]  ^  matrix  destroying  [ME]  in  the 
process.  The  Gaussian  Elimination  Method  with  column  pivoting 
is  used  in  determination  of  the  inverse  of  the  matrix.  Since 
this  subroutine  is  small  and  performs  only  the  one  function, 
it  has  not  been  flow  charted  but  is  listed  in  Appendix  C. 

C.  PROGRAM  LIMITATIONS 

Infinite  summations  such  as  those  which  appear  in  equa¬ 
tions  (36)  and  (41)  ,  can  not  be  performed  on  a  digital  com¬ 
puter  which  deals  only  with  finite  numbers.  Therefore  infinity 
must  be  replaced  by  a  finite  approximation.  In  the  case  of 
equation  (36),  an  analysis  of  G]j£x(an)  I  versus  n  was  con¬ 
ducted  with  results  shown  in  Figure  6.  As  shown  in  this  figure, 
when  n  reaches  a  value  of  about  10  the  expression  is  down  by 

-3 

at  least  a  factor  of  10  relative  to  its  largest  value  close 
to  the  origin.  The  value  n  =  20  is  sufficiently  large  to 
be  considered  infinite  for  all  practical  purposes.  Therefore 
equation  (36)  becomes 


00 


l 

n=0 


GJEI2 


G.  (0 )  20 

-V-*  l.  Van'S)IW 

n=i 


(45) 
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Figure  6.  GilEx(an)|2  vs  n  with  er  =  2.2  and 

D/A  =  0.02  for  Five  Arbitrary  Values 
of  A '/A 
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After  establishing  this  limit  on  the  summation,  computer 
errors  (overflows)  occurred  at  various  points  in  the  summa¬ 
tion.  An  investigation  revealed  that  the  problem  originated 
in  the  dyadic  Green's  function  component  when  became 

too  large.  Varying  the  size  of  h^  produced  overflows  only 
when  h  was  increased  beyond  a  certain  point;  a  situation 
where  the  walls  around  the  central  slotline-like  structure 
became  far  removed  from  the  field.  Limits  were  able  to  be 
placed  on  y^h^  such  that  if  its  value  became  greater  than  some 
value  X  in  this  case  X  »  80,  then  y^h^  was  set  equal  to  X 
because  of  the  limited  interaction  between  the  walls  and  the 
field.  When  the  field-wall  interaction  becomes  a  factor, 
the  overflow  problem  no  longer  exists.  Limitations  of  this 
kind  are  found  in  the  computations  of  the  wavelength  as  well 
as  the  impedance. 

The  high  valued  arguments  of  the  hyperbolic  functions  just 

discussed  cause  additional  problems  when  determining  the 
0  ll 

coefficients  A  through  D  in  the  power  flow  computations. 

In  the  inversion  of  the  matrix  [ME] ,  these  large  values  become 
small  ones  and  when  squared  produce  underflows  in  the  computer 

system  since  the  magnitude  of  the  elements  were  less  than 

•78  0 

10  .  Limitations  were  placed  only  on  the  coefficients  A  , 

0  h 

D  and  D  which  caused  the  underflow  problem.  The  coefficients 
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were  limited  to  magnitudes  of  5*10  .  This  problem  was 

compounded  by  the  fact  that  the  coefficients  are  complex  and 
they  may  also  take  on  either  plus  or  minus  signs. 


During  the  investigation  of  the  overflow  problem,  it 
was  noted  that  when  W/b  =1,  E  =  0  for  n  >  0  and  so  G, (0) 

X  x 

was  the  only  important  term  in  the  computation.  Therefore 

whenever  W/b  =  1  only  the  n  =  0  term  is  computed  in  both 

the  wavelength  and  impedance  computations. 

In  the  characteristic  impedance  equation  (41) ,  the 

coefficients  in  this  infinite  series  were  found  to  decay 

rapidly  so  that  a  finite  approximation  yields  good  results. 

Figure  7  shows  the  characteristic  impedance  as  a  function  of 

N,  the  number  of  terms  in  the  truncated  series,  for  various 

slot  widths.  Termination  of  the  impedance  computation  occurs 

when  the  difference  between  two  consecutive  values  of  Z  in 

o 

the  iteration  loop  falls  below  0.1%  or  when  the  upper  limit 
of  N  *  50  is  reached.  A  minimum  of  seven  iterations  is  set 
for  the  larger  values  of  W/b  except  for  the  situation  when 
W/b  =1.  In  this  case,  only  the  N  =  0  term  is  non-zero. 

D.  INPUT/OUTPUT  DATA 

Proper  input  data  formatting  is  essential  for  any  success¬ 
ful  computer  program.  Input  data,  similar  to  that  shown  in 
Tables  I  and  II  are  required  with  the  program  'FIN-LINE'.  A 
minimum  of  three  data  cards  is  necessary  to  make  up  the 
input  data  deck.  The  first  card  specifies  the  dimensions  of 
the  shielding  and  the  permittivity  of  the  dielectric  material 
(b/D,  hj/D,  l^/D  and  er) .  This  information  is  entered  in  a 
floating  point  format  with  the  column  requirements  specified 
in  the  program  listing  in  Appendix  C.  Note  that  for  a  ridged 


Characteristic  Impedance,  Z  ,  vs.  N,  the 
number  of  terms  in  the  truncated  series, 
with  e  =  2.2  and  D/X  =  0.0212  for  a  WR(19) 
shield1 


13.8 

13.3 

17.8 

1.0 

2 

0.1 

0.0170 

15 

0.  0 

0.5 

0. 0212 

24 

0.5 

0. 0021 
0.0021 


OOOGOOOOO  0  00000000000000000  OflDGOOOGO  0  OOGOOO 
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Table  I.  Sample  Input  Data  for  a  Ridged  Waveguide 
Structure  with  WR(19)  Shield 


19.8  18.3  17.8  2.2 

3 

0.02  0.0170  15  0.0  0.  0021 

°*5  0-0170  15  0.0  0.04 

1-0  0.01  15  0.0  0.01 


Table  II.  Sample  Input  Data  for  Fin-line  Structure 
with  WR (19 )  Shield 
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waveguide,  the  dielectric  substrate  thickness  is  zero.  This 
would  result  in  infinite  values  for  the  first  three  values 
of  the  input  data  (b/D,  h^/D  and  h^/D) •  To  eliminate  this 
problem  select  any  thickness  greater  than  zero,  but  let  its 
dielectric  constant  equal  one. 

The  second  card  specifies  the  number  of  initial  value 
data  cards  that  follow  with  nine  being  the  maximum  allowed. 
Sample  input  of  Table  I  indicates  that  two  cards  will  follow 
the  second  card  and  sample  input  of  Table  II  shows  three  data 
cards  follow. 

The  third  card  needed  for  the  input  contains  initial 
values  for  the  computations  that  follow.  Six  pieces  of  informa¬ 
tion  are  supplied  by  this  card.  The  first  two  pieces  of 
information  contain  the  initial  values  of  W/b  and  D/A .  Next 
a  two  digit  number  specifies  the  number  of  different  values 
of  W/b  (first  digit)  and  D/A  (second  digit)  that  will  be  com¬ 
puted.  On  the  last  card  of  Table  I,  the  number  24  indicates 
that  the  wavelength  and  impedance  will  be  determined  for  two 
different  values  of  W/b  and  for  four  different  values  of  D/A 
for  each  W/b.  The  last  two  values  are  the  incremental  values 
of  W/b  and  D/A,  respectively.  The  value  of  0.5  will  be  used 
to  increment  the  initial  value  of  W/b  so  that  the  second  value 
will  be  W/b  *  1.  The  value  0.0021  will  be  used  to  increment 
the  initial  value  of  D/A  so  that  the  four  values  of  D/A  are 
0.0212,  0.0233,  0.0254  and  0.0275. 

The  data  shown  in  Tables  III  and  IV  is  the  computer  output 
from  execution  of  the  program  listed  in  Appendix  C  with  the 
input  data  of  Table  I  and  Table  II,  respectively. 


SLOT -TRANSMISSION  LIN'S  ON  A  DIELECTRIC  SU8STRATC 
WITH  SHIELDEC  WALLS 
(FIN-LINE) 


DIMENSIONS:  6/0*18.80  Hl/0*18.80  H2/D=17.80  EPSR*  1.0 


RIDGED  WAVEGUIOE  STRUCTURE  EXISTS 
FOR  THE  FOLLOWING  VALUES  OF  w/B 


W/B*  0.100  0/L*  0.0170  LP/L*  1.1559  ZO*  215.7  OHMS 
W/B*  0.100  0/L*  0.0191  LP/L*  1.1176  Z0=  208.5  OHMS 
W/B*  0.100  C/L*  0.0212  LP/L*  1.0918  ZO*  203.8  OHMS 
W/B*  C.100  0/L*  0.0233  LP/L*  1.0738  ZO*  200.5  OHMS 
W/B*  0.100  0/L*  0.0254  LP/L*  1.0613  ZO*  198.1  OHMS 
W/B*  0.500  C/L*  0.0212  LP/L*  1.1988  ZO*  382.6  OHMS 
W/B*  0.500  D/L*  0.0233  LP/L*  1.1559  ZO*  369.1  OHMS 
W/B*  0.500  C/L*  0.0254  LP/L*  1.1262  ZO*  359.6  OHMS 


TE-20  EMPTY  GUIDE  MODS  PROPAGATES  FOR 
VALUES  OF  0/L  GREATER  THAN  0.027 


W/3*  0.500  D/L*  0.0275  LP/L*  1.1051  ZO*  352.5  OHMS 

RECTANGULAR  WAVEGUIOE  STRUCTURE 
EXISTS  FOR  W/B*  1.0 


W/B*  1.000  0/L*  0.0212  LP/L*  1.2840  ZO*  484.1  OHMS 
W/B*  1.000  0/L*  0.0233  LP/L*  1.2176  ZO*  459.2  OHMS 
W/B*  1.000  0/L*  0.0254  LP/L*  1.1738  ZO*  442.4  OHMS 
W/B*  1.000  D/L*  0.0275  LP/L*  1.1426  ZO*  430.6  OHMS 


Table  III.  Output  Data  Obtained  Using  Sample 
Input  Data  from  Table  I 
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slot-transhissijnhlin|ponca  dielectric  substrate 

(FIN-LINE) 


DIMENSIONS:  0/0=18.80 


Hl/D-18.80 


H2/0= 1 7. 80 


SPSR*  2.2 


1=  0.020  C/L  =  0.0170 
3*  0.020  D/L=  0.0191 
3  =  0.020  0/L*  0.0212 
)=  0.020  0/L*  0.0233 
3  =  0.020  0/L=  0.0254 


LP/L=  0.9183 
LP/L=  0.9035 
LP/L=  0.8933 
LP/L=  0.8855 
LP/L=  0.8793 


3=  0.500  C/L=  0.0170  LP/L=  1.2457  Z0= 

TE-20  EMPTY  GUIDE  MODE  PROPAGATES  FOP 
VALUES  OF  D/L  GREATER  THAN  0.027 

}  =  0.500  D/L*  0.0570  LP/L*  0.9465  Z0= 

3  =  0.500  D/L*  0.0970  LP/L*  0.9043  Z0= 

3=  0.500  D/L*  0.1370  LP/L*  0.8683  Z0= 

SURFACE  WAVES  PROPAGATE  FOR  VALUES  OF  ! 
GREATER  THAN  0.169;  RUM  TERMINATED 

DIELECTRIC  SLAB-LOADED  WAVEGUIDE 
STRUCTURE  EXISTS  FOR  W/8*  1.0 

CUTOFF  FREQUENCY  OF  THE  DOMINANT  MODE 
OF  EMPTY  GUIDE  IS  C/L=  0.013 


138.6  OHMS 

138.1  OHMS 

138.2  OHMS 
138.9  OHMS 
139.8  OHMS 

421.2  OHMS 


501.2  OHMS 
700.0  OHMS 
776.9  OHMS 


1.000 


OF  EMPTY  GUIDE  IS  C/L*  0.013 
C/L*  0.0100  LP/L*  ****** 


3=  1.000  D/L*  0.0200  LP/L*  1.2660 


1.000  0/L*  0.0300 
1.000  0/L*  0.0400 


LP/L*  1.0707 
LP/L*  1.0191 


(*  1.000  C/L*  0.0500  LP/L*  0.9957 


***** 

508.3  OHMS 
464.9  OHMS 
493.8  OHMS 
553.6  OHMS 


Table  IV.  Output  Data  Obtained  Using  Sample 
Sample  Input  Data  from  Table  II 
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In  addition  to  the  tabulated  results  of  the  computations 
for  the  wavelength  and  impedance,  message  statements  are  also 
printed  under  certain  conditions.  All  the  messages  that  can 
be  generated  by  the  program  are  found  in  the  output  of  Tables 
III  and  IV. 

Since  various  structures  may  take  form  by  varying  the 
input  parameters,  description  of  the  type  of  structure  that 
exists  will  be  printed.  For  ridged  waveguide,  the  permittivity, 
er»  must  be  one  and  W/b  <  1  whereas  the  structure  becomes 
empty  rectangular  waveguide  when  W/b  =  1  (shown  in  Table  III) . 
With  er  >  1  and  W/b  <  1,  a  fin-line  structure  is  formed.  Since 
this  analysis  is  of  fin-line  structure,  no  output  statement 
is  printed.  When  W/b  =  1  and  e  >  1,  a  dielectric  slab  loaded 
waveguide  structure  exists  (Table  IV) . 

The  analysis  of  the  fin-line  structure  presented  here 
concerns  itself  mainly  with  operation  in  the  dominant  TE1Q 
limit  mode.  At  the  lower  end  of  the  dominant  mode  frequency 
range  when  D/X  <  D/X^,  the  cutoff  frequency  of  the  empty  wave¬ 
guide,  an  output  statement  with  this  information  is  printed. 

The  cutoff  wavelength  of  the  empty  waveguide,  X^,  is  equal  to 
2a  where  a  =  h^+l^+D.  Since  the  structure  takes  other  forms 
than  just  an  empty  waveguide,  operation  can  exist  below  D/Xc. 
With  the  presence  of  a  dielectric  slab  in  the  fin-line  and 
slab  loaded  waveguide  operation  exists  below  D/xc.  Therefore 
analysis  below  the  cutoff  frequency  of  empty  waveguide  is 
possible  but  a  lower  limit  is  set.  When  the  value  of  x'/x 
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exceeds  6,  indicating  the  approach  of  cutoff,  computations 
terminate  without  printing  values  of  \ ' /\  and  ZQ  for  the 
values  of  D/x  for  which  x '/x  >  6. 

At  the  high  end  of  the  dominant  mode  frequency  range  when 
X  =  a,  the  TE20  of  the  empty  guide  can  propagate.  Therefore 
a  statement  is  printed  when  D/x  >_  D/a.  For  a  value  of 
D/X  >_  l/(4/T),  surface  waves  propagate  on  the  metal-backed 
dielectric  slab  and  normal  operation  is  no  longer  possible 
so  the  run  is  terminated  for  that  value  of  D/x. 

The  total  program  length  including  all  subroutines  is 
48  K  bytes  of  storage.  The  sample  computer  runs  presented 
here  with  outputs  shown  in  Tables  III  and  IV  took  4  and  5 
minutes,  respectively  which  includes  compilation  time  of 
about  20  to  30  seconds. 
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IV.  NUMERICAL  RESULTS  AND  COMPARISONS 


To  check  the  accuracy  of  the  numerical  results  generated 
by  the  computer  program,  comparisons  were  made  with  data 
available  in  the  literature  for  various  fin-line  sub-structures. 
These  structures  vary  from  ridged  waveguide  at  one  end  to 
dielectric  slab  loaded  waveguide  at  the  other. 

A.  RIDGED  WAVEGUIDE 

As  mentioned  earlier,  the  fin-line  structure  shown  in 
Figure  1  becomes  ridged  waveguide  with  zero  thickness  ridges 
when  =  1  and  W/b  <  1  or  when  the  dielectric  substrate 
thickness  D  is  reduced  to  zero . 

The  design  of  ridged  waveguide  was  first  treated  by  Cohn 
[Ref.  20]  and  later  by  Hopfer  [Ref.  21]  and  Lager lof  [Ref. 

22] .  The  wavelength  and  impedance  of  the  ridged  waveguide 
structure  can  be  given  by  the  following  expressions 


AV 

A 


[1  -  (-r-)  2]1/2 


Z 
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[1 


Z 

O00 


-  (t^)2]1/2 


(46a) 


(46b) 


The  cutoff  wavelength  Ac  was  obtained  from  the  curves  pre¬ 
sented  in  [Ref.  21]  for  several  values  of  W/b.  Using  this 
information  and  equation  (46a),  the  wavelength  ratio  A'/A  was 
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determined.  These  calculated  values  were  then  compared  with 
the  results  using  the  spectral  domain  method  of  the  computer 
program  discussed  in  Section  III.  Figure  8  shows  the  two 
results  with  agreement  to  within  0.5%  where  the  discrete 
points  are  from  equation  (46a)  . 

A  similar  comparison  of  the  impedance  was  performed  using 
equation  (46b)  where  ZQ  was  obtained  from  curves  found  in 
[Ref.  22]  and  the  numerical  results  of  the  program  'FIN-LINE! 
The  results  appear  in  Figure  9  where  again  very  good  agreement 
exists  with  a  difference  of  only  2%. 

B.  DIELECTRIC  SLAB-LOADED  WAVEGUIDE 

A  variation  of  the  fin-line  structure  of  Figure  1  where 
W/b  -  1  and  er  >  1  results  in  a  dielectric  slab  loaded 
rectangular  waveguide.  This  structure  has  been  studied  in 
considerable  detail  by  Vartanian,  et  al. ,  [Ref.  23].  In  their 
work,  the  dielectric  slab  was  centered  in  the  waveguide  and 
they  presented  an  analytical  expression  for  the  (voltage) 
impedance,  Z^v  (using  the  notation  of  [23])  at  the  center  of 
the  slab.  The  impedance  computed  in  this  analysis  was  speci¬ 
fied  at  the  edge  of  the  dielectric  slab.  Therefore  a  relation¬ 
ship  between  the  two  impedances  must  be  defined  before  any 
comparison  can  be  made.  The  impedance  at  the  edge  of  the  slab 
can  be  easily  defined  by  the  expression 


68 


eT  *  1.0 


where  is  the  field  at  the  edge  of  the  slab  and  E^ntr 

is  the  field  at  the  center  of  the  slab.  Thus, 

=  Z  cos 2 { (47b) 
o  pv  2  s 

where  the  various  quantities  are  defined  in  [22]  as 

q2  =  (2tt)2  (^->2  [  e  -  (X/X')2] 

X  r 

s  =  c/2 

and  C  =  D  is  the  dielectric  slab  thickness.  The  characteristic 
impedance  of  a  slab  loaded  waveguide  was  computed  using 
equation  (47b)  (where  Z^v  was  obtained  from  Vartanian's 
curves)  and  the  spectral  domain  method.  The  results  of  the 
computations  are  compared  in  Figure  10  where  the  discrete 
points  are  calculated  from  equation  (47b) .  It  can  be  seen 
that  the  results  agree  very  closely.  The  difference  is  less 
than  1%. 

The  wavelength  of  the  slab  loaded  waveguide  may  be 
determined  analytically  by  the  transverse  resonance  proce¬ 
dure  [Ref.  16]  which  is  found  by  writing  the  field  expressions 
for  an  exact  solution  of  the  problem,  applying  boundary  condi¬ 
tions  at  the  air-dielectric  interfaces  and  thereby  obtaining 
the  determinantal  equation  for  the  structure.  The  wavelength 
ratio  is  found  by  determining  the  value  of  x'/x  which  satisfies 
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the  determinantal  equation.  Figure  11  shows  the  variation 
of  the  wavelength  ratio  with  dielectric  thickness,  D,  for 
several  values  of  er  obtained  by  the  spectral  domain  method. 

The  edge  of  the  slab  was  placed  at  the  mid-point  of  the  broad 
wall  of  a  WR(19)  waveguide  operating  at  40  GHz.  The  dis¬ 
crete  points  are  those  values  determined  by  the  transverse 
resonance  approach.  Excellent  agreement  was  obtained  between 
the  results  of  the  two  methods.  The  greatest  difference  between 
the  different  approaches  was  only  0.6%. 

C .  SLOTLINE 

If  W/D  <_  2  and  is  sufficiently  high  for  the  fin-line 
structure  of  Figure  1,  the  presence  of  the  shield  will  have 
little  effect  if  the  walls  are  sufficiently  far  removed  from 
the  slot.  In  this  case  the  structure  will  behave  like  a  slot¬ 
line.  This  behavior  of  fin-line  is  illustrated  in  Figure  12 
where  the  wavelength  and  characteristic  impedance  is  plotted 
for  er  *  20  and  a  slot  width  normalized  to  the  dielectric 
thickness,  W/D  =  1.  The  discrete  points  represent  the  wave¬ 
length  (  □  )  and  the  impedance  (  )  of  a  slotline  obtained 

from  [Ref.  24]  for  the  same  parameters  er  and  W/D  as  the  fin- 
line.  Agreement  between  the  fin-line  and  the  slotline  is 
within  1%  for  the  wavelength  and  5%  for  the  impedance.  The 
impedance  curve  for  a  ridged  waveguide  structure  with  the 
same  fin  separation  is  plotted  for  reference. 
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NORMALIZED  FREQUENCY  D/X 


Figure  12.  Wavelength  Ratio,  and  Characteristic 

Impedance,  Z  vs.  Normalized  Frequency  D/A 
for  Fin-line°and  Slotline.  Fin-line  Has 
WR (19)  Shield  and  Fins  Centered 


D.  MEASUREMENTS 


Very  little  experimental  data  is  available  in  the 
literature  for  any  significant  comparison.  However , fin-line 
wavelength  measurements  have  been  presented  by  Meier  [Ref.  1] . 
The  guide  wavelength  in  the  millimeter  fin-line  structure  was 
measured  across  the  band  of  26.5  -  40  GHz  by  a  sliding  short- 
circuit  technique.  The  fin-line  used  a  D  =  10  mil  substrate 
with  er  ■  2.2.  A  WR (28)  shield  (a  =  .280",  b  =  .140")  was 
used  along  with  a  fin  separation  W/b  =  .128  (W  =  .018") . 

The  results  of  this  experimental  work  are  plotted  in  Figure 
13.  A  computer  run  was  made  for  a  structure  with  the 
parameters  reported  by  Meier  and  the  results  were  compared 
with  Meier's  measured  values.  This  comparison  is  illustrated 
in  Figure  13  where  an  agreement  to  within  1%  can  be  seen. 


Figure  13.  Computed  and  Measured  Values  of  Wavelength 
Ratio  X'/X  vs.  Frequency  for  Fin-line 
with  WR(28)  Shield.  Fins  are  Centered 


V.  FIN-LINE  DESIGN  CURVES 


A  problem  exists  with  the  preparation  of  design  curves 
for  structures  like  fin-line  because  of  the  number  of  inde¬ 
pendently  variable  parameters  which  describe  the  structure. 

This  difficulty  can  be  alleviated  to  some  degree  by  practi¬ 
cal  considerations.  First,  fin-lines  are  constructed  with 
a  shield  that  is  compatible  with  the  dimensions  of  standard 
rectangular  waveguide  for  the  millimeter  wavebands.  Above 
22  GHz,  all  the  standard  waveguide  structures  have  aspect 
ratios,  b/a  =  0.5.  Furthermore,  the  fins  are  normally 
printed  on  5  mil  thick  Rogers  Duroid  5880  dielectric  substrates 
with  er  =  2.2.  Design  curves  have  been  provided,  which 
appear  in  Figures  14-17,  for  structures  with  parameters  dis¬ 
cussed  above  for  the  26.5  -  40  GHz,  40  -  60  GHz,  60  -  90  GHz 
and  the  90  -  140  GHz  waveguide  bands. 

It  can  be  seen  from  Figures  14-17  that  even  with  fin 
separations  of  a  few  mils,  low  values  of  impedance  are  diffi¬ 
cult  to  achieve.  In  some  applications  lower  impedances  than 
those  obtained  here  which  are  on  the  order  of  125  -  150  ohms, 
are  desirable.  One  method  of  accomplishing  this  objective 
would  be  to  increase  the  dielectric  constant.  This  can  be 
seen  by  comparing  Figure  12  for  a  W/D  =  1  and  Figure  15 
for  W/b  =  .053.  For  this  W  =  5  mil  slot,  it  is  found  that 
the  impedance  decreases  from  about  165  ohms  to  90  ohms  when 
the  dielectric  constant  is  increased  from  er  =  2.2  to  er  =  20. 
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Wavelength  Ratio  X '/A  and  Characteristic 
Impedance  Z  vs.  Frequency  for  Fin-line  with 
WR(28)  Shield.  Fins  are  Centered 
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Wavelength  Ratio  X'/x  and  Characteristic 
Impedance  Z  vs.  Frequency  for  Fin-line  with 
WR (12)  Shield.  Fins  are  Centered 
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The  disadvantage  of  this  method  is,  of  course,  a  reduction 
in  the  wavelength  by  a  factor  of  2.5  causing  the  already 
small  circuit  dimension  to  become  even  smaller. 

Another  possible  way  of  achieving  a  lower  value  of  impe¬ 
dance  may  be  in  the  relocation  of  the  fins  toward  the  side- 
walls  of  the  guide.  Figure  18  illustrates  the  wavelength 
and  the  characteristic  impedance  of  a  fin-line  with  the  fins 
shifted  midway  between  the  center  and  the  sidewall  of  a  WR(19) 
waveguide.  A  comparison  of  Figures  18  and  15  indicates  a 
2-4%  change  in  the  wavelength  of  the  guide  for  W/b  =  1  and 
less  for  smaller  values  of  W/b. 

The  values  of  the  characteristic  impedance  for  small 
values  of  W/b  are  relatively  unchanged  since  the  fin-line 
structure  behaves  like  slotline  and  the  presence  of  the  shield 
has  little  affect.  For  values  of  W/b  >  0.5,  noticeable 
decreases  in  the  impedance  have  taken  place.  The  most  signi¬ 
ficant  decrease  occurs  for  W/b  =  1  where  the  impedance 
decreases  from  500  ohms  for  the  centered  dielectric  slab  to 
about  275  ohms  for  the  off-centered  configuration.  In  this 
condition  the  structure  is  a  slab  loaded  waveguide  and  since 
the  low  dielectric  constant  results  in  little  change  in  the 
empty  guide  fields,  the  results  are  as  expected  based  on  the 
voltage  impedance  definition  for  the  empty  rectangular  wave¬ 
guide.  A  point  of  interest  is  the  impedance  for  W/b  =  1 
falling  below  that  for  W/b  =  0.5.  This  indicates  that  for 
the  off-centered  configuration,  the  impedance  increases  with 
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FREQUENCY  (GHZ) 


Figure  18.  Wavelength  Ratio  \'/X  and  Characteristic 
Impedance  Z  vs.  Frequency  for  Fin-line 
with  WR(19)  Shield.  Fins  are  Located 
Half  Way  between  Side  Wall  and  a  Guide 
Center 


W/b  until  a  maximum  is  reached  and  then  decreases  until 
W/b  =  1.  Overall,  the  relocation  of  the  fins  toward  one  side 
of  the  waveguide  appears  not  to  have  achieved  the  lowering 
of  the  impedance  for  small  values  of  W/b  as  anticipated. 
Therefore,  of  the  two  methods  discussed,  increasing  the 
dielectric  constant  accomplishes  the  lowering  of  the  char¬ 
acteristic  impedance  desired;  but  as  mentioned  before,  the 
circuit  dimension  shrinkage  is  unavoidably  linked  to  this 
action. 
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VI.  CONCLUSIONS 

A  spectral  domain  technique  for  analysis  of  the  dominant 
mode  of  fin-line  has  been  presented.  It  has  been  shown  that 
this  efficient  numerical  method  can  be  implemented  compu¬ 
tationally  to  accurately  calculate  the  dispersion  character¬ 
istic  and  characteristic  impedance  as  a  function  of  frequency. 
A  matrix  formulation  of  the  problem  permitted  the  elements 
of  the  dyadic  Green's  function  to  be  computed  and  circumvented 
the  extensive  algebraic  manipulations  associated  with  the 
formulation  of  the  resulting  eauations.  Numerical  result? 
obtained  by  using  this  method  have  been  presented  and  compared 
with  existing  data  for  various  structures.  Excellent  agree¬ 
ment  was  obtained  in  all  cases  thus  establishing  accuracy  of 
the  method  and  the  versatility  of  the  computer  program. 

Design  curves  for  millimeter-wave  fin-lines  of  practical 
interest  were  also  included  here.  Both  the  centered  and 
off-centered  fin  configurations  were  discussed  and  the  off- 
centered  fin  location  was  shown  to  result  in  no  significant 
impedance  change  for  small  values  of  W/b. 

In  this  thesis,  particular  interest  was  devoted  to  the 
analysis  of  fin-line.  The  method  employed,  however,  is  com¬ 
pletely  versatile  since  it  is  a  general  formulation  of  an 
inhomogeneous  transmission  media  problem.  Therefore,  it  can 
be  clearly  seen  that  by  varying  the  parameters  of  the  fin-line 
structure,  it  may  exhibit  the  characteristics  of  ridged 
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waveguide,  dielectric  slab  loaded  waveguide,  slotline,  and 
conventional  rectangular  waveguide.  All  of  these  structures 
are  fin-line  substructures  and  as  such  are  included  in  the 
approach  presented. 


87 


APPENDIX  A 


Derivation  of  Spectral  Domain  Matrices 
[Me]  and  [Mj] 

The  continuity  conditions  given  in  equations  (9)  through 
(18)  are  transformed  via  (24)  into  the  two  dimensional 
Fourier  domain  which  yields 
at  v  =  D 


kci 


kc2  *!«vd) 


e,  3*1  e  3»2 

*  -jc.nr»2(an,D)-j„U5-^(„n,D) 


kcl  *1<VD) 


Et <<0 

z  n 


e  3$i 

-j«nr$i(an'D)-i^37-(“n'D)  = 


E  (a  ) 
x  n 


(A-l) 


kcl‘l(“n'D>-kc242(“n'D)  =  W 


h  3«T  h 

-janr$1(oln,D)+jwei— (an,D)+janr$2(an,D) 


3  <I>2 

•jue23y-(an'D)  =  Jz(an) 


at  y  =  0 


*e2(*n,o) 
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e  .  0  *2 

■j ^n1" $ 2  ( an' ^  “j - (an-0)  = 


;?h 

_:lV}3(an'0’_:)“,,,3F(V0) 


k=2  ‘X'0) 


kc3  ’^n'0’ 


(A-2 ) 


h  3$2  h  3$T 

■janr$2(an'0,+j“e2^r(an'0)  =  *j  anF  *3  (an '  0)  +: j  “e3  W~{  V 0) 


where  y=y^=y2=u3. 

The  solutions  to  the  two  Helmholtz  equations  given  by 
equations  (26a) — ( 26f )  are  now  substituted;  the  dependence 

Q  h 

of  the  coefficients  A  through  D  on  a  is  understood  and 
from  here  on  is  not  shown  explicitly. 


(k^sinhy-j^hj^)  Ae- (k^2sinhy2D)  Be- (k^2coshy2D)  Ce  =  0 


(A—  3 ) 


-  ( janrsinhy1h1)  A®+  ( j  c^Tsinhy^)  Be+  ( j  c^TCOShy^)  Ce+  ( j  wy  y^sinhy^h^) A 


+  ( j  aip  y2COShy2D)  B^+ ( j  wy  y 2sinhy2D)  C^1  =  0 


(kc2)Ce_(kc3  sinhy3h2)De  =  0 


(janr)Ce+(janrsinhy3h2)De-(ja)yy2)Bh 


+  ( jwy y3sinhy3h2) D  = 


<ko2>ch 


(k^3coshy3h2) Dh  =  0 


(A-4) 


(A-5) 


(A-6) 


(A-7) 


(A-8) 


( j  U)£2y2^  ^  j  we  2V  3COS^y3^2^  ^e~  ( 3 <*nr )  C*1 


■  ^  T  T  ^*2 


+  ( janrC0Shy3h2) Dr 


=  o 


-  ( j  a  rsinhYi  h, )  Ae+ (i  am  y-,  sinhy-.  h. )  A^1  =  E  (a  ) 
n  11  1  11  xn 


(A-9 ) 


(k^sinhY^h^)  Ae  =  Ez(<*n) 


(A-10 ) 


(k^coshY-j^^)  Ah-(k^2sinhY2D)  Bh_  ^kc2COshv2D)  ch  =  ^x^n'1 


(A-ll) 


0 

-  ( jaje1Y1coshY1h1)  A  -  ( j(i)£2Y2coshY2D)  Bc-  ( j  we2Y2sinhY2D)  C 


e  ,  . 


(A-12) 


-(ja  rcoshY-,h,  )  Ah+(  ja  rsinhY->D)  Bh+  ( j  a  rcoshY-,D)Ch  =  J  (a  ) 
n  11  n  2  n  2  zn 


This  set  of  linear  equations  (A- 3)  through  (A-12)  may 
be  written  in  matrix  form,  such  that 


where  the  8*8  matrix  is  [M^]  of  equation  (29a),  obtained 
by  combining  (A-3)  through  (A-10) . 

The  combination  of  equations  (A-3)  through  (A-8) ,  (A-ll) 
and  (A-12)  form  the  matrix  equation  (29b)  repeated  below  in 
expanded  form. 

mll  m12  m13  0  0  0  0  0 

m21  m22  m23  0  m25  m26  m27  0  B 

0  0  m33  m34  0  °  0  °  Ce 

0  0  m.-  m,,  0  m,-  0  m.Q  De 

43  44  46  48 

0  0  0  0  0  0  m57  m58  Ah 

0  m62  0  m64  0  0  m67  m68  B" 

0  0  0  0  m^5  m^g  m^?  0  Ch 

m81  m82  m83  0  m85  m86  m87  0 

The  matrix  elements  of  [M^]  and  [Mj]  are  normalized  at 
this  point  with  respect  to  D,  the  dielectric  substrate 
thickness;  but  before  describing  them,  the  normalized  con¬ 
stants  used  in  these  elements  are  presented. 

(YlD)2  =  (Y3D)2  =  (ctnD)  2-K2tt)  2  [  (  A/A  ’  )  2  -  1]  (D/X)2 

(Y2D)2  =  UnD)2  +  (2tt)  2  [  (A/A  * )  2  -  er]  (D/A)2 
2 

If  (y^D)  <  0  then  the  following  is  defined 
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(  y1”D)  2  =  -(y-jD)2  =  -(y3D)2 
(y2"D)2  =  -(y2D)2 
Other  normalized  constants  are 

(kclD)2  =  (kc3D)2  =  (2fl)2f1-(x/x')2nD/X)2 
(kc2D)2  =  (  2tt)  2  [  er-  (  \/X  * )  2  ]  (D/X)2 

uyD  a  240  tt 2  (D/X) 


oie  =  oie  3D  =  1/60 (D/X) 

ijj  €  2  D  =  er/60(D/X) 


0D  =  2 it  (D/X)  (X/X  ' ) 


anD 


n2ir  (D/b)  $  even 

(2n-l)  it  (D/b)  <f>h  odd. 


The  normalized  matrix  elements  are  now  presented  in  two 

2 

forms.  The  first  of  the  element  equations  is  for  (y^D)  >  0 

2 

and  the  second  is  for  (y^D)  <  0.  For  the  matrices  [M^l 


and  [Mj], 


the  elements  through  m68  are  the  same. 
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m25 

m26 

m27 

m28 

m31 

m33 

m34 

m35 

m41 


j  UuD)  (vj^DJsinht  U-jD)  (hj/D)  ] 

■ 

-j  UuD)  (Yl"D)  sint  (Yl"D)  (h1/D)  ] 

* 

j  (uigD)  (^DJcosh^D 
- (uuD)  (y2"D) cosy  2 " D 
j UuD) (y^D) sinhy^D 

< 

-j UuD) (y2"D) sinY2"D 

0 


m32  =  0 


(kc2D) 


(kc2D) 


- (kc3D)  sinht (y3D) <h2/D) ] 


-j(kc3D)2sin[(Y3,,D)  (h2/D) 


m36  =  m37  =  m38  =  0 


m42  =  0 


(y3d)2  1  0 
(y3d)2  <  0 
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m4  3 


It\  a  A 

44 


UnD)  (3D) 

(anD) (SD) 

-  ( a n_D )  (  8D)  sinhf  (y  3D)  (hj/D)  ] 
-j (anD) (8D)sin[(y3"D) (h2/D) ] 


m45  =  0 


m46 


m4  7 


-j  (ojuD)  (Y2D) 
(uuD)  ( y 2 " D) 


m48 


m51 


m57 


j (upD)  (y3D) sinh[ (y3D)  (h2/D)  ] 

-j (wpD) (y3"D) sin[ (y3"D) (h2/D) 3 
=  m52  =  m53  =  m54  =  m55  =  m5g  =  0 

,kc2D)2 

,  <kc2D»2 
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F/t  1/3 


-(kc3D)2cosht (y3D) (h2/D) ] 
-(kc3D)2cos[(Y3HD) (h2/D)3 


m71  =  < 


(a^D)  (  SD)  sinhf  (y^)  (h^D)  ] 


j  UnD)  ( SD)  sin  [  (  y1"D)  (hx/D)  ] 


E  _  E  E 
m72  m73  -  m74  =  0 


m 


75 


j  UyD)  (y.j^D)  sinh  [  iy^D)  (h^/D)  ] 


-j  UuD)  (Yl"D)  sin[  (y1"D)  (hj/D)  ] 


m76  m77  m78  _  0 


m 


81 


(kclD)  2sinh  [  (y^D)  (hj/D)  ] 


j(kclD)":sin[(Y1"D)  (hj/D)] 


=  mf_  =  m_.  -  = 


m82  m83  m84  m85 


m71  m72  m73  -  m74 
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0 


m87  m88 


=  0 


m75  - 


(kclD) ^cosh [ ( Y]D) (h1/D) ] 


(kclD)'!cos[(Y1HD)  (h1/D)  ] 


2 

~(kc2D)  sinh^D 

"j  <kc2D) 2sinY2"D 

2 

-(kc2D)  coshy2D 
- (kC2°) 2cosy2"D 

-j (uie^D) (y^D) cosh[ (ylD) (h^/D) ] 
(uie^D)  (yj/'D)  cos  [  (y-^’D)  (h^/D)  ] 
-j  (uie2D)  (Y2D)coshY2D 
(u^D)  {y2"D)  COSY2"D 
-j (uejD) (y2°) sinhY2D 
j  (o)e2D)  (y2"D)  sinY2"D 

(anD)  (6D)cosh[  (y^D)  (h^/D)  ] 

(anD)  (8D)  cos  [  (y1"D)  (hj/D)  ] 


J 

m86 


m 


87 


m 


88 


(a  D) ( 8D) cosh (  (y,  D) (h./D) ] 
n  ■».  J. 

. 

(ctnD)  ( 8D)  cos  [  (  y1"D)  (hx/D)  ] 

- (a„D) ( 8D) sinhy-D 
n  i 

-j (anD) (8D)siny2HD 


Once  the  matirx  elements  are  determined  for  the  existing 

condition  (y^  *  0) ,  the  [Mg]  1  is  created  for  the  purpose  of 

0 

eliminating  the  unknown  coefficients  A  through  D  in  (31) . 
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APPENDIX  B 


Time  Average  Power  Flow  in  the  Spectral 
Domain  of  Fin-line 


The  following  expressions  for  the  time  average  power  flow 
in  the  three  spatial  regions  are  derived  from  equation  (40) . 
The  solutions  to  Helmholtz's  equations  given  in  (26a)  through 
(26f)  are  substituted  in  the  transformed  scalar  potentials 
of  (41) .  In  the  following  derivations  the  dependence  of  the 

0  h 

coefficients  A  through  D  on  the  variable  a  is  understood. 

n 

Also  to  be  noted  is  that  y^  =  Y3  since  regions  1  and  3  are 
assumed  air  filled. 

Region  1 

2 

Case  la:  (^D)  >_  0,  Yl  real 


=  Ae  sinh  Y-j^D+h^y) 


=  Ah  cosh  Y^D+h^-y) 


°»1  e 

J—  =  -y-jA  cosh  Y1(D+h1-y) 


-Y.jAh  sinh  y1(D+h1~y) 


The  term  Pla  will  be  used  for  the  power  flow  in  region  1  for 
2 

the  case  (y^)  >  0.  Substituting  the  Helmholtz  equation 

solutions  for  region  1,  (40)  becomes 


i 
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D+hl 

Pla  =  JS  Re  E  /  {  t-8a)eia^|Ae|2-6^u1Yj!Ah!2+jS2cln>1Ae(Ah)  * 

n— — D 

-jkJanYi(Ae)*Ah]  [sinh2y1(£Hh1-y)  ] 

+  [-ea)M1a2|Ah|2-3a)e1Y2|Ae|2-t-j62anY1Ae(Ah)  * 

“jkl0lnYlAh(Ae)*l  tcoah^tD+hj-y)  ]}dy  .  (B-l) 

Carrying  out  the  integration 
D+h1 

/  sinh2Y1 (D+h^-y) dy  =  j^-tsinh  Uy^)  ^Yjh^ 

J  cosh2Y1(D+h1-y)dy  =  ^-[sinh  ( 2Y]h1) +2Y1h1] 


Equation  (B-l)  becomes 


00  0u)  e .  a 


Pla  "  15  Re  l  {C - Ti~I1|Ae|2-e<,U1Y1|Ah|2+j62anAe(Ah)* 

n— — ®  T 1  1 


-  jkJan(Ae)*Ah] tsinh(2Y1h1)-2Ylh1] 

2 

h  2  „  i  „e ,  2 .  _2  .e,.h^ 

+  — - - 1 A  |  -6ue^Y1|A  I  +36  aRA  (A  )* 


■jk2anAh(Ae) *}  tsinh(2Y1h1)+2Ylh1] } 


(B-2) 


101 


Equation  (B-2)  is  the  final  result  for  the  power  flow 

2 

xn  region  1  when  (y^D)  1.  0.  Putting  (B-2)  in  normalized 

forms  produces  the  following  result. 


'  la 


i  n  "  (6D  a„Dr  .e  , 

-  ifH)  Re  Y  i  f _ 1  n  i  A  1 2 

8  b  Re  1  11  (YlD) 

n=-«  1  '  D 


Ah  i  2 


— jl  +(3D)  (wu^D)  (y-^D)  |— j 


-j(6D)2(a  D)^(-^-*+j  (k,D)2(ct  D)^- 


h,  h,  ( 3D)  Uu.D)  (a  D)  h  - 

l  n  I A  1 1 


[sinh2(YlD)  (-^)-2(YlD)  (-^)]  +  [- 


Ty^dT 


+  (3D)  (ue.D)  ( Yi  D)  ] 2- j  ( SD)  2  (a  D)  £  •(^)- 
1  1  D2  n  D2  D2 


+j  (k1D)2(anD)^  [sinh2(YlD)  (^)+2(Y;lD)  (^)  ]  }  (B-3) 


2 

Case  lb:  (y^D)  <  0,  yi  is  imaginary  in  which  case 

( Y 2_ ” D )  2  =  -  (YiD)  2 

=  jAe  sin  Yj^'^D+hj^-y) 

♦  ^  *  A*1  cos  Yl"(D+h^-y) 

3*® 

__  .  -  j  Y^"Ae  cos  Y1"(D+h1-y) 
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«  *h 


ajT  =  Y x "A  sin  Y1" (D+h1-y) 

Substituting  the  values  into  equation  (41) ,  integrating  and 
normalizing  with  respect  to  D  results  in 


lb 


Ah  i  2 


l  n  “  (8D)  (we,  D)  (a  D)  .e  ,  .. 

-  |(g)  Re  l  {[ - (y-  „by  n - |^|2+(BD)(u>u1D)(Yi"D)|^- 

n=-»  T 1  D  D 


-(BD)2(a  D)^  (k1D)2(a  D)-^- 

n  D2  D2  1  n  D2  D2 


h.  h.  (BD)  (wu,  D)  ( a  D) A  -h  , 

[2(Yl"D)(^)-sin2(Yl"D)(1l)]  +  [ - {y4b)  |A  '2 


6  0  h  * 

+  (SD)  (we-jD)  (Yi"D)  l^-|2+(BD)2(anD)^j  ^ 


+  (k1D)2(anD)-^-  ^][2(Yl"D)  (^)+sin2(Yl"D)  (^)  ]  }  (B-4) 


This  is  the  final  result  for  the  power  flow  in  region  1 
2 

when  (YlD)  <  0. 

A  shorter  version  of  the  derivation  will  be  used  for  the 
remaining  cases  wherein  the  regional  cases  will  be  described 
and  the  final  normalized  result  will  be  shown.  This  will  be 
similar  to  that  done  for  case  lb. 


For  region  2,  the  power  flow  expressions  are  the  same  as 
for  unshielded  slowline  except  that  the  Fouriar  integral  is 


replaced  by  a  summation  and  the  interval  2n  is  replaced  by 
the  interval  b. 

2 

Case  2a:  (y2D)  1  y2  is  real 

6  6  6 
$2  =  B  sinh  Y2y  +  C  cosh  Y2y 

<J>2  *  Bh  sinh  Y2y  +  Ch  cosh  Y2y 

=  y 2  [Be  cosh  Y2y  +  c  sinh  Y2y] 

3 

°  U.  U 

•Jy= ■  =  y2[B  cosh  Y2y  +  cn  sinh  Y2y] 

Substituting  into  equation  (41)  and  integrating  followed  by 
normalizing  with  respect  to  D  produces  the  following. 


!,£,  V  fr<BD)^«2D)(anD)  |Be|2  ( 3D)  ( uu  2D)  ( ^D) 

8  b  n-—  {y2D)  V  +  ^2D) 


Ce  i  2 . 


+  (SD)  («c2D)  (y2D)  |-j|  + ( SD)  (uu2D)  (y2D)  |=j 


Ch  |  2 


+j  IlDl’l.  D]  (5|  (k  D)  2  („  D) 

D  D  D2  2  n 

.  (JbV  C^+Bh  (C®)* 

D2  D2  ?  “^2~ 


„  .  (SD) (ue,D) (anD)2  - 

•  [sinh(2Y2D)-2y2D)]+[ - - 2 - |^|2 


(SD) (uw,D) (a  D) 2  _h  0 
_ 2  n  I  c  i  2 


(Y2D) 


Cn  f  2  pe  o 

~j[  +(SD)  ( to e _D)  (y„D)  |^*-| 

D  •  ^  DZ 


*  (8D)(u„2D)(y ,D)|^|2+j(SD,2(a  d)(1B^c!+b!  _(CV 

D  "  D2  D2 

-j  (k2D)  (otnD)  (-J  -^2 —  +  “^2 —  [sinh(2y2D)+2y2D] 

‘“I”  *»’”’>  IjJ  ^ 

+j2(BD)2(a  D)  (5|  IcV).j2(k  d)2(o  d, 

D  i  D^  D2  2  n 

r*^  /o®  \  *  ti h 


./(Cc)  C, (Be)  B 

V——*—  — *T - x -  — 


D  D  D  D 


-y] tcosh (2y2D) -1] } 


This  is  the  final  result  for  the  power  flow  in  region  2  when 


Case  2b: 


2 

(y2D)  <  0,  Y2  is  imaginary  in  which  case 

(y2"D) 2  =  -(y2D52 

<J>2  =  jBe  sin ( Y2"y)  +  Ce  cos(Y2"y) 

$2  -  jBh  sin(Y2"y)  +  Ch  cos(Y2"y) 


3<tr 

=  jY2"B  cos  Y2"y  -  y2"  c  sin  Y2"y 

3  4*5  h  h 

=  jY2"B  cos  Y2Hy  -  Y2"cn  sinY2"y 

2 

The  final  result  for  region  2  when  (y2D)  <  0  after  integration 

and  normalization  follows. 
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ao 


2b 


v  UP)  U^D)  CnD)  Be  ,«DIU!D)|lWi.1k 

8<b'  J-.  <v2"D)  ?  + - (T75T -  ^ 


p®  9  r  * 

+  (6D)  (ue2D)  (y2’’D)  |=yP+(SD)  (u^D)  (y2”D) 


Ch  2 


+  (3D)2(V»(^  IcV.iiV^ 

D  D  D 


+  (k2D)2(anD)  (-^J-  [2Y2"D-sin(2Y2"D)] 


D  D“  D“  D* 
2 


(6D) (we-D) (a  D)  e  ,  (8D) (wu,D)(a „D)  _h  _ 
1  ( y2"D)  1  2  '  +  '  ■  ,,T"  I  *  1 


(y2"D) 


i  Be i 2  . 


Bh  i  2 


+  (8D)  (we2D)  (y2"D)  |— jj  +(6D)  (uju2D)  (y2"D)  |— y. 


+  (8D)2(a  D)  (  -By  -^  )  +  (k  D)  2  ( a  D)  (5-  i£-l - (iL  1  £_)  j 

D  D  D*  D*  z  n  dz  £)Z  jj*  g* 

(a  D)2-(y  ’’D)2 

•  [2Y2"D+sin(2Y2"D)  ]+j  [  (3D)  (ue2D)  ( — ,lT^  2 - ) 


.Be  _(Ce)_(Be)  C* 

\  n  n  ^  ‘  “J  + 


2  2 
D  D  D‘ 


(y2"d) 


^nD)2-(Y2"D)2 


r~  ^2)+(6D)  (“^d)  <— Ty2-d) — ) 
h  IcV_(Bh'  *  ~h 

D?  D2  D"  D 

2 


(!l  TT-TT-  5l)+2(BD)2(anD) 


e  *  li  6  *  h 

+  2  (k2D)  ‘  (a  D)  (-^J-  ]  [1-cos  (2Y'’D)  ]  } 

n  D  D  D2  D2  2 


(B-6) 
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Region  3 

Case  3a:  (y3D)2  _>  0,  y3  is  real 

$3  =  D  sinh  y3(h2+y) 

$3  »  Dh  cosh  y3(h2+y) 

a$3  e 

Ty  =  Y3De  C°Sh  Y3(h2+y) 

a<fth 

a  ^  3  h 

Jy~  =  y3D  sinh  Y3(h2+y) 

Substituting  these  expressions  into  (41)  ,  integrating  and 
normalizing  yields 


3a 


ID  *  (BD)  (<oe,D)  (<*nD>*  _e  , 

-£<£)  Re  l  {{— - 1 - - - '  D  i  2 


n=-« 


(y3D) 


ily|‘'+(SD)  (uu  ,D)  (y-,D)  l^yl2 
D  J  D 


+  j(8D)2(a  D)  ^^-j{k3D)2(a  D)iiLi 
n  D*2  3  n  d2  D2 

,  .  h2  h,  <6D)  <»ii-D)  (c.  D)  2  _h  , 

•tsinh2(v3D)  <-£>-2<y  D)  <-£>1  +  1 - p-3-D-)--B - 1^1 

w3  '  D 

+  (8D)  (o»e,D)  (y.D)  |5y|2+j  (<?D)2(a  D)2»  1*U- 

D  n  d1  DZ 


2  (Del  *  n^1  h_  h_ 

-3(k3D)  (anD)-i-^i-  Hy]  Esinh2(Y3D)  {-^)+2{y3D>  (^-)  )  } 


h„ 

D 


(B-7) 
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This  is  the  final  result  for  region  3  power  flow  for 
(y3D)2  ^  0. 

2 

Case  3b:  (y3D)  <  0,  y3  is  imaginary  in  which  case 

( Y  3 " D ) 2  =  -(y3D)2. 


=  jDe  sin  y3"(h2+y) 


$3  =  Dh  cos  y3"(h2+y) 


=  3y3"D  cos  y3" (h2+y) 


-Y3"Dh  sin  y3" (h2+y) 


Substituting  into  equation  (41),  integrating  from  y  =  ~h2 
to  y  =  0  and  normalizing  the  result  with  respect  to  D  produces 
the  following. 


ID  r  <BD)  (“nD)  ne  9  nh  9 

Re  l  H - - - - l^l  +(6D)  (o,w3D)  (y  ”D)  |^.|2 

n=-°°  T3  D  J  D 


+  (3D)2(a  D)^  i°^l_+(k  D)2(a  D)i2j_  p\ 
n  u  u  J  n  v  u 

h2  h-  (6D)  (o)m,D)  (a  D) 2  h  , 

•[2(y  »D)(-|)-sin2(y3"D)  (-^)]+[ - r-^,nr  5 - |2y|2 


(y3"d) 


+  (8D)  (ue.D)  (y  "D)  \£\ 2-  (  SD)  2  (a  D)  ^ 

D  n  D*  D4 

1  (n®\  *  h-  h- 

-(k3D)  (anD)-^-  H_]  [2(y3"D)  (-^+sin2(y3"D)  (-^)  J } 


(B-8) 


This  is  the  final  normalized  expression  for  the  power  flow 

2 

in  region  3  under  the  (y^D)  <  0  condition.  It  may  be 

noted  that  a  minus  sign  appears  in  front  of  all  the  equations, 
this  is  because  a  wave  traveling  in  the  -z  direction  was 
assumed. 
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FINE-LINE 


APPENDIX  C 


Computer  Program  'FIN-LINE' 
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